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Annual  Report  -  AFOSR  Grant  83-0367 
Research  Period:  1  Oct.  1983  -  30  Sept.  1984 

A.T.  Mattick  and  A.  Hertzberg 


Abstract 


Technical  results  for  the  first  year  of  research  on  AFOSR  grant  83-0367 
are  summarized.  This  work  involved  a  theoretical  study  of  radiation  trans¬ 
fer  in  droplet  clouds,  for  predicting  the  performance  of  the  liquid  droplet 
radiator.  The  effects  on  droplet  sheet  emissivity  of  non-isotropic 
scattering  by  droplets  and  of  the  development  of  a  temperature  profile  on 
the  droplet  sheet  are  shown  to  be  small. 


-1- 


1.  Introduction 

This  report  summarizes  the  results  of  a  theoretical  analysis  of  radia¬ 
tion  transfer  in  droplet  clouds  carried  out  during  the  period  1  Oct.  1983  to 
30  Sept.  1984  on  AFOSR  grant  83-0367.  This  program  of  research  was  ini¬ 
tiated  to  establish  a  theoretical  foundation  for  predicting  the  performance 
of  the  liquid  droplet  radiator  (LDR),  an  advanced  heat  rejecttion  system  for 
space  power.  A  schematic  of  the  baseline  configuration  for  this  radiation 
is  shown  in  Fig.  1.  Waste  heat  is  transferred  to  a  liquid,  which  is  sprayed 
into  space  in  the  form  of  submillimeter  droplets.  The  droplet  trajectories 
must  be  controlled  accurately  enough  that  all  the  liquid  can  be  collected, 
once  the  droplets  have  cooled  by  radiating  their  energy  to  space.  This  con¬ 
cept  has  the  promise  of  being  much  lighter  than  conventional  radiator 
designs  which  use  solid  surfaces  to  reject  waste  heat. 

Previous  analyses  of  the  radiation  by  the  rectangular  sheet  of  droplets 
carried  out  by  the  investigators^1’^) ,  had  assumed  that  the  temperature  did 
not  vary  in  the  sheet  thickness  direction  and  that  the  droplets  scattered 
radiation  isotropically.  The  chief  aims  of  the  present  program  were  to 
derive  a  formalism  for  radiative  transfer  in  an  LDR  droplet  sheet  which  in¬ 
cluded  the  effects  of  non-isotropic  scattering  and  a  time  dependent  tempera¬ 
ture  profile,  and  to  find  numerical  solutions  for  the  droplet  sheet 
emissivity.  The  investigators  have  accomplished  these  goals  during  the 
first  year  of  the  program,  and  have  initiated  an  experimental  program  to 
test  the  predictions  of  the  theory,  by  measuring  the  radiation  emitted  by  a 
cloud  of  droplets.  air  fo-t  c.r"'rrz  nr  ''.tf'M'tftc  (Ajsc) 
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The  theoretical  analysis  of  radiation  transfer  including  non-isotropic 
scattering  indicated  that  for  droplet  sheet  parameters  of  interest  for  the 
LDR,  non-isotropic  and  isotropic  results  do  not  differ  substantially.  Like¬ 
wise  allowance  for  a  non-uniform  temperature  profile  in  the  sheet  thickness 
direction  has  only  a  small  effect  on  the  power  emitted  by  an  LDR  droplet 
sheet. 

Investigation  of  a  solar  thermophotovoltaic  system  for  space  power  was 
terminated  because  an  engineering  analysis  of  the  liquid  silicon  storage 
scheme  for  this  system  indicated  that  materials  incompatibility  problems  are 
insurmountable.  Another  research  area  of  the  proposal  for  this  grant  was 
the  investigation  of  the  emissive  properties  of  metals  and  alloys,  along 
with  investigation  of  methods  for  increasing  emissivities  of  low  emissivity 
materials  which  otherwise  have  promise  as  candidate  LDR  fluids.  We  carried 
out  a  literature  survey  of  the  emissive  properties  of  metals,  both  to  find 
optical  characteristics  of  classes  of  metals  and  to  find  a  suitable  theory 
for  predicting  emissivities  of  metal  mixtures.  This  extensive  survey 
revealed  that  the  existing  data  on  metal  emissivities  (particularly  for  tin 
and  tin  alloys,  and  for  alkali  metals)  was  very  scarce,  particularly  for  the 
liquid  phase,  and  secondly,  that  there  is  at  present  no  theory  of  metal 
emissivity  which  can  acccurately  predict  the  temperature  dependence  of 
emissivity.  We  concluded  that  further  efforts  in  this  area,  either  in 
measuring  liquid  metal  emissivities,  or  in  making  fundamental  contributions 
in  the  theory  of  optical  properties  of  metals,  was  beyond  the  scope  of  the 
program. 
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2.  Analysis  of  Radiation  Transfer  in  Droplet  Clouds 

The  purpose  of  the  theoretical  program  was  to  develop  a  more  complete 
model  of  radiation  transfer  in  droplet  clouds  used  for  the  liquid  droplet 
radiator  (LDR)  than  had  previously  been  used.  Our  prior  work  on  this 
problem  had  assumed  a  simple  model  of  an  LDR  droplet  sheet  in  which  the 
droplet  temperature  was  uniform  and  in  which  scattering  of  radiation  by 
droplets  was  isotropic.1*^  The  radiation  interaction  with  the  droplets  was 
treated  using  geometrical  optics.  The  principal  modifications  to  this  model 
addressed  in  the  present  program  were  the  allowance  for  a  non-uniform 
temperature  profile  and  for  non-isotropic  scattering  by  the  droplets.  Our 
goal  was  to  determine  the  effects  these  phenomena  would  have  on  the 
radiative  characteristics  of  a  droplet  sheet  as  compared  to  the  results  from 
our  simple  model . 

Model  of  a  Planar  LDR  Droplet  Sheet 

In  the  baseline  configuration  of  the  liquid  droplet  radiator  (Fig.  1.) 
a  rectangular  array  of  heated  liquid  droplets  is  ejected  into  space  at  the 
droplet  generator.  The  width  W  of  this  array  (y-di rection)  is  generally 
much  larger  than  the  thickness  S  (x  direction):  The  droplets  traverse  a 
distance  L  (sheet  length,  z  direction),  radiating  waste  heat  into  a 
radiation  vacuum  (no  incident  radiation). 

The  droplet  sheet  is  analogous  to  a  solid  radiating  surface  of  a 
conventional  radiator,  and  it  is  conveniently  characterized  by  a 
hemispherical  emissivity  g$  from  its  "surfaces,"  which  are  the  top  and 
bottom  of  a  rectangular  sheet  and  have  dimensions  W*L.  Generally  W  and  L 
are  of  comparable  magnitude  and  the  area,  WL,  is  chosen  to  meet  the  heat 
rejection  requirements  of  a  specific  application,  i.e.,  ^rad  =  ^es°B^ 


where  og  is  the  Stefan-Boltzmann  constant.  The  sheet  emissivity  es  may  vary 
with  distance  z  along  the  direction  of  droplet  propagation.  The  sheet  emis¬ 
sivity  depends  on  the  number  density  of  droplets  and  sheet  thickness,  which 
are  assumed  constant  along  the  sheet  length,  and  also  on  the  absorption  and 
scattering  properties  of  individual  droplets  which  may  in  turn  depend  on 
temperature.  The  variation  in  sheet  emissivity  thus  arises  from  a  decrease 
in  the  mean  temperature  of  the  droplets  (in  an  x-y  slice  through  the  sheet) 
as  they  traverse  the  transit  length  L,  and  also  from  the  development  of  a 
temperature  profile  through  the  thickness  of  the  sheet  (x-di rection) . 
Although  the  droplets  have  a  uniform  temperature  as  they  are  ejected  into 
space,  a  temperature  variation  in  the  thickness  direction  will  develop 
because  the  droplets  nearer  the  surfaces  will  cool  faster  than  interior 
droplets . 

Our  analytical  model  of  radiation  transfer  in  a  droplet  sheet  makes  two 
simplifying  assumptions  based  on  the  relative  sheet  dimensions.  The  sheet 
thickness  S  is  chosen  so  that  the  optical  depth  tqx  =  naS  of  the  sheet  in 
the  thickness  direction  will  be  of  order  unity,  (n  is  the  number  density  of 

p 

droplets  and  a  =  na^  is  the  frontal  area  of  a  droplet.)  For  practical 
designs  (see  Ref  3  )  this  dimension  will  be  on  the  order  of  centimeters, 
whereas  the  dimensions  W  and  L  will  be  on  the  order  of  10's  of  meters.  Edge 
effects  will  be  important  only  over  distances  corresponding  to  ~1  optical 
depth.  Thus  variation  of  temperature  in  the  y-direction  (sheet  width)  will 
only  be  significant  over  a  fraction  ~0.001  of  the  sheet  width.  Our  first 
assumption,  then,  is  that  the  sheet  is  effectively  infinitely  wide,  and 
there  is  no  variation  of  droplet  temperature  or  radiation  intensity  in  the 
y-direction.  A  second  simplification  arises  from  the  fact  that  S/L  <<<  1 
and  that  for  a  practical  LDR  design  the  fractional  temperature  decrease  of 
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droplets  in  flight  is  <0.5.  The  logarithmic  derivative  of  temperature  with 
optical  depth  along  the  z  direction  is 

1  J]_  ~  aT/T  =  aT/T  S 
^  3tz  toz  tox  ^ 

and  is  of  order  10"^. 

Droplets  at  a  given  point  in  the  sheet  only  interact  radiatively  with 
those  within  an  optical  distance  of  order  unity.  Since  the  fractional 
variation  in  droplet  temperature  (and  droplet  radiation)  over  one  optical 
depth  in  the  z  direction  is  extremely  small,  we  neglect  it.  The  second 
assumption,  then,  is  that  the  temperature  is  uniform  in  the  z-direction. 

The  two  assumptions  reduce  the  model  to  a  single  dimension,  x,  and  the 
problem  is  to  determine  the  temperature  evolution  of  an  infinite  slab  of 
radiating  droplets,  which  initially  have  a  uniform  temperature.  This 
infinite  slab  represents  a  slice  of  the  droplet  sheet  as  it  propagates  from 
the  droplet  generator  to  the  collector.  The  above  considerations  assure 
that  this  model  will  predict  the  radiative  characteristics  of  an  LDR  droplet 
sheet  (finite  W  and  L)  very  precisely. 

Radiative  Transfer  Equations 

In  the  model  described  above,  the  infinite  droplet  sheet  will  have  a 
uniform  initial  temperature,  but  since  the  droplets  nearer  the  surface  of 
the  sheet  cool  more  rapidly  than  interior  droplets  a  non-uniform  temperature 
profile  will  evolve  normal  to  the  sheet,  i.e.,  T  =  T(x,t).  It  is 
nonetheless  useful  to  characterize  the  sheet  by  a  hemispherical  sheet  emis- 
sivity  based  on  the  average  sheet  temperature: 


T(t)  =  i  /  dx  T(x,t) 
^  o 


(1) 
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The  sheet  emissivity  is  defined  to  be: 
es(t)  =  J(t)/0BT4(t) 


(2) 


where  J ( t )  is  the  hemispherical  power/unit  area  emitted  by  a  surface  of  the 
sheet.  This  convention  is  useful  because  the  energy  balance  equation  then 
has  the  same  form  as  that  for  a  uniform  temperature  system: 


m  r  dT  - 
msC  dt  “ 


2esaBT4 


(3) 


where  m$  is  the  fluid  mass/ (projected  sheet  area),  and  C  is  the  specific 
heat  of  the  fluid.  The  principal  aims  of  our  program  were  to  determine  the 
time  evolution  of  e$  and  to  evaluate  the  effect  of  non-isotropic  scattering 
on  es.  When  considering  fluids  which  exhibit  structure  in  their  absorption 
and  emission  spectra,  and  also  for  the  spectral  variation  in  droplet 
scattering  due  to  wave  effects,  it  is  often  necessary  to  consider  the 
spectral  sheet  emissivity 

eSv(t)  = 

where  Jv  is  the  sheet  hemispherical  intensity/(unit  frequency),  and  J® ( t )  is 
the  blackbody  hemispherical  i ntensity/(uni t  frequency): 

Jy(T)  =  -1  v3/(exp(hv/kl)-l)  . 
c c 


The  relationship  between  eSv  and  eg  is  given  by: 
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es  =  /dveSvJ®/aBT4  (4) 

and,  of  course, 

J(t)  =  /dvJv(t) 
aBT4  =  /dvJ®(T)  . 


The  droplet  temperature  is  governed  by  the  equation  of  radiative 
transfer: 


al  (x,n,t) 

- - - cose  =  no  {-  Iv(x,n,t)(eov+«ov) 


+  e0vJ5[T(x,t)]/w 


+  ^/dn’Iv(x,n'  ,t)ov(n'+ft)} 


and  by  radiative  energy  balance  for  a  droplet: 


mC  HUiil.  /dv  {-4o  covJ?[T(x,t)] 


+  e0va/da  Iu(x,fi,t)}  . 


(5) 


(6) 


In  these  equations,  n  is  the  number  density  of  droplets,  m  is  the  droplet 
mass,  C  is  the  fluid  specific  heat,  o  =  tt is  the  projected  geometrical 
cross-section  of  a  droplet  (a  =  droplet  radius),  eQv  is  the  frequency- 
dependent  droplet  emittance  (often  called  the  absorption  efficiency)  and 
is  the  scattering  efficiency.  Note  that  the  total  scattering  and  absorption 
cross-sections  at  frequency  v  are  ouv  and  oe0v,  respectively.  The  spectral 
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intensity  per  solid  angle,  I  (x,Si,t),  in  droplet  sheet  is  not  only  a 
function  of  frequency,  position  and  time,  but  also  of  direction  Q  =  (e,<^) 
where  0  is  the  polar  angle  from  the  sheet  normal  x,  and  i|»  is  the  azimuthal 
angle.  In  fact,  Iv  is  independent  of  t Jj  due  to  the  symmetry  of  the  problem. 
The  non-isotropy  of  scattering  is  expressed  by  the  phase  function  ov  which 
is  normalized  to  4tt ,  i.e.,  /da' $(&'-*■&)  =  4tt .  (For  isotropic  scattering 
$v  =  1.)  The  quantities  eQv,  u>v,  and  <t>v  can  be  obtained  from  the  Mie 
scattering  theory  if  the  complex  index  of  refraction,  n  =  nj+in£,  is  known 
as  a  function  of  frequency.  These  quantities  are  generally  complicated 
functions  of  the  size  parameter  z  =  2ira/X.  For  droplets  which  are  large 
compared  to  the  radiation  wavelength  (z  >100)  the  scattering  exhibits  very 
rapid  oscillations  as  a  function  of  scattering  angle  and  a  very  strong  peak 
over  a  narrow  range  of  angles  in  the  forward  direction.  However,  if 
averaged  over  a  small  range  of  frequencies  (or  equivalently,  z),  the 
scattering  is  very  nearly  isotropic,  aside  from  the  forward  peak.  This  peak 
can  be  subtracted  from  the  scattering  distribution  because,  for  the  present 
problem,  forward  scattered  light  does  not  differ  from  transmitted  light. 

For  smaller  values  of  z,  the  non-isotropy  varies  more  slowly  with  angle  and 
cannot  be  "averaged  out."  In  fact,  an  average  of  $  over  a  range  of 
wavelengths  large  enough  to  incorporate  2  or  more  nodes,  still  yields  a  non¬ 
isotropic  distribution.  This  fact  was  the  primary  basis  for  our 
investigation  of  the  effects  of  non-isotropic  scattering  on  droplet  sheet 
emi ssi vity. 

Owing  to  the  symmetry  of  the  problem  about  a  pane  midway  through  the 
sheet  (at  x  =  S/2),  we  have  ( I v(x,0>^,t)  =  Iv(S-x,n-0 ,t) .  The  formalism  is 
also  simplified  if  we  change  from  independent  variable  x  to  a  new  variable 
tv  defined  below,  and  introduce  the  relative  efficiencies  eQv,  uQv: 
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tv  =  nax(e0v+u)v)  =  T(e0v+uv) 


e0v  ~  e0v^e0v+“v^ 


w  =  u)  /  i  +to  )  . 

v  v'  '  Ov  v ' 


For  convenience  the  tilde  will  be  omitted  below.  We  can  formally  integrate 
eq.  (5) ,  to  obtain: 


MV-*)  ■  ^  6XP[' 

0 


Tv“Tv  i  .  iB 


]  •  {eovJv^T(Tv»t)]/Ti  (8) 


+  uj  de 'sine 1  [lv(t^,e 1 ,t)Y+(e,e' )+iv(Tov-r;,e' ,t)Y_(e,e' )]} 


where 


2  7T  ov(e.0;e',*) 

Y±(e,e‘)  =  ^  /  d*  { 

0  «j>v(e,0;n-e',4-) 


In  order  to  isolate  the  effects  of  a  nonuniform  temperature  profile  and 
non-isotropic  scattering  on  droplet  sheet  emissivity,  simplified  approaches 
will  be  taken  -  isotropic  scattering  was  assumed  in  determining  the 
evolution  of  a  temperature  profile,  and  a  uniform  temperature  was  assumed  in 
determining  the  effects  of  non-isotropic  scattering.  We  shall  describe  the 
analysis  of  each  case  separately. 

Temperature  Profile  of  a  Droplet  Sheet 

In  the  determination  of  the  temperature  profile  of  a  droplet  sheet,  we 
assume  isotropic  scattering,  so  that  #  =  1  and  y-=y+=1/2.  We  further 

assume  that  the  droplets  are  gray  bodies,  that  is,  the  scattering  and 
absorption  efficiencies  have  no  frequency  dependence.  With  these 


•*.  . iv.  w--.  ^ '"o 


assumptions  we  can  immediately  integrate  eq.  (8)  over  frequency,  and  carry 
out  the  frequency  integral  in  eq.  (6)  to  obtain: 


I(x,p,t)  =  -  /d-r'expf-  J  .  {e0ogT4(r '  ,t)/n 

p  0  u 

1 

+  j  /dp' [I (x '  ,p' ,t)  +  I (xQ-x ' ,p 1 ,t)]j  (9) 

o 

mC  =  -4e0aaBT4(T,t) 

ot 

1 

+  2tt  e0a  /dp[I(x,p,t)  +  I(xQ-x,p,t)]  (10) 

0 

where  u  =  cose.  Denoting  the  quantity  in  curly  braces  by  f(x‘,t), 
introducing  a  dimensionless  variable  f(x,t)  =  ? (x ,t ) •  (Tr/e0agl4) ,  where  TQ  is 
the  initial  temperature,  and  also  using  a  reduced  temperature  f  =  T/TQ,  we 
obtai n : 


-  t. 


f(x,t)  =  T  (x,t)  +  -y  j  dx'  E^(  t-x'  )f(x‘,t) 


(ID 


3-^-V-£)-  =  -T4(?,t)  +  eQf  (x  ,t) 
3t  u 


(12) 


where  the  time  variable  has  been  replaced  by  the  dimensionless  quantity: 


t  =  t  4e0aogTQ/mC<]j 


(13) 


and  we  have  reincorporated  the  tildes  for  the  quantities  x,  w,  and  eQ. 

co 

En(x)  is  the  exponential  integral  /  dz  e"xz/zn. 
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The  coupled  equations  11  and  12  are  to  be  solved  to  obtain  the  tempera¬ 
ture  profile  f(x,t),  with  the  initial  condition  T(x,0)  =  1.  The  solution 
parameters  are  the  modified  optical  depth  xQ  =  x0(eQ+2)  and  the  modified 
droplet  emittance  lQ  -  e0/(e0+u),  since  to  =  l-eQ.  For  large  droplets 
e0+to  =  1,  so  that  x=x  and  e0  =  eQ.  Assuming  this  to  be  the  case,  we  can 
again  drop  the  tilde  notation.  The  hemispherical  emissivity  of  the  droplet 
sheet  is: 

1 

e$(t)  =  /udu  I(x0,y,t)/aBT4(t) 
o 

or 

e$  =  /°dx  E2(x)f(x,t)/T4(t)  .  (14) 

n  o 

Equations  11  and  12  were  solved  in  an  iterative  fashion.  The  intensity 
function  f  was  determined  for  a  given  distribution  T  by  solving  integral 
equation  11  numerically.  The  approach  was  similar  to  that  used  in  our 
earlier,  uniform  temperature  calculations  where  the  half  sheet  (since  the 
problem  is  symmetrical)  was  broken  up  into  individual  stations.  The 
integral  equation  was  then  transformed  into  a  matrix  equation: 


fi(t)  =  f|(t)  +  |  lAijfj(t)  (15) 

j 

and  f  could  be  determined  by  a  simple  matrix  inversion.  The  distribution 
f .,•  ( t )  was  then  used  as  a  source  function  to  advance  f^(t)  by  a  small  time 
increment  at  using  eq.  (12). 


Numerical  solutions  for  the  temperature  profile  as  a  function  of  time 
were  obtained  for  a  wide  range  of  optical  depths  and  droplet  emittances. 

From  an  initial  state  of  uniform  temperature  there  is  a  period  of  rapid 
adjustment  of  the  temperature  profile  whereby  the  temperature  of  droplets 
near  the  surface  of  the  sheet  falls  more  rapidly  than  the  interior  tempera¬ 
ture.  For  optical  depths  between  0.5  and  4  and  for  emittances  between  0.1 
and  1  this  adjustment  was  nearly  complete  at  reduced  time  t~l.  After  this 
adjustment  period,  although  the  absolute  temperature  decreased  steadily  with 
time,  the  relative  temperature  distribution  T(x)/T(x=S/2)  remained  nearly 
constant.  This  distribution  was  dependent,  of  course,  on  the  optical  depth 
and  droplet  emittance.  The  larger  the  optical  depth  the  greater  the 
fractional  temperature  decrease  from  the  center  to  the  edge  of  the  sheet. 

Figure  2  illustrates  the  magnitude  of  the  temperature  variations  normal 
to  the  sheet.  Plotted  are  curves  showing  the  temperature  (normalized  to  the 
temperature  at  the  center  of  the  sheet)  as  a  function  of  optical  distance 
from  the  center  of  the  sheet,  for  combinations  of  (e0,t0)  of  (0.8,1), 

(0.2, 2. 5),  and  (0.1,5),  and  for  a  time  t*  at  which  the  mean  temperature  has 
decreased  to  0.8  of  its  iniital  value  T0.  These  parameters  are  typical  of 
those  which  might  be  used  in  a  liquid  droplet  radiator,  and  the  profiles 
represent  the  temperature  distribution  at  the  droplet  collector  for  an  LDR 
designed  so  that  the  mean  droplet  temperature  decreases  by  20%  during 
flight.  For  eQ  =  0.1,  xQ  =  5,  the  surface  temperature  is  12%  below  the 
center  temperature,  and  for  the  other  cases  the  difference  is  6-7%. 

Although  these  differences  are  substantial,  considering  that  the  fractional 
decrease  in  mean  temperature  is  20%,  the  effect  on  sheet  emissivity  is  not 
large.  This  is  illustrated  in  Fig.  3,  which  shows  the  emissivity  of  the 
sheet  as  a  function  of  time,  normalized  to  t*.  For  an  LDR  these  curves 


-13- 


correspond  to  sheet  emissivity  as  a  function  of  distance  from  the  droplet 
generator  to  the  droplet  collector.  For  eQ=0.1,  tq=5,  the  emissivity  falls 
from  0.486  to  0.466,  a  fractional  decrease  of  4%,  while  for  eo=0.8,  tq=1, 
the  fractional  decrease  in  es  is  only  ~2%.  The  assumption  of  a  uniform  tem¬ 
perature  profile  would,  of  course,  yield  a  constant  emissivity  equal  to  the 


initial  emissivity  of  the  non-uniform  temperature  case. 

We  conclude  that  for  values  of  eQ  and  tq  of  interest  for  LDR 
applications,  it  is  reasonably  accurate  to  assume  a  constant  emissivity  of 
the  droplet  sheet,  computed  as  though  no  temperature  profile  develops.  Our 
earlier  investigations'*  have  shown  that  a  reasonable  compromise  between 
minimizing  the  sheet  area  (larger  optical  depth)  and  minimizing  biioet  mass 
(smaller  optical  depth)  would  result  from  choosing  tq  =  0.5/eQ.  This  choice 
results  in  a  sheet  emissivity  (es  =  0.5±0.05  for  droplet  emittances  eQ 
betwen  0.05  and  1.  The  existence  of  a  temperature  profile  may  of  course  be 
critically  important  when  the  heat  transfer  medium  operates  close  to  the 
freezing  point,  in  which  case  the  droplets  near  the  sheet  surface  may 
solidify.  On  the  other  hand,  in  cases  where  the  fractional  decrease  in  mean 
droplet  temperature  is  small,  a  temperature  profile  may  not  have  time  to 
fully  develop. 

Effects  of  Non-Isotropic  Scattering 


As  mentioned  previously,  to  isolate  the  effects  of  non-isotropic 
scattering  on  droplet  sheet  emissivity,  computations  of  this  emissivity  were 
made  assuming  a  uniform  temperature  droplet  sheet.  Thus  the  source  function 
J®(T)  i°  ^e  integrand  of  the  integral  equation  for  I  (eq.  (8))  is  time- 
independent,  as  is  Iv  itself.  Denoting  the  function  in  curly  braces  in 
eq.  (8)  by  Lv  and  defining  a  new  dimensionless  intensity  function 
Lv  =  Lv/e0j|J(T),  we  obtain: 


*  C  UV  V  V  <•  O".'  vV.\  ■  *  '  '  •  • 


i.  i’ V?7"VV”iT?Wly???7??lTiT?crCVWi,viSVs?r»W?i;rrri  i%?fntff?fr<r(rCr(rkfCf(rCrtr(rCB^^P»wWr#Sfvi 


-14- 


Lv(Tv»e) 


ir/2 

1  +  co/  de'sine1 
o 


1 

cose 1 


(16) 


+ 


{  Y*(8,e'}  /  it;  Lv (t^, , e 1  )exp[- 
o 


Tv"Tv 

cose* 


] 


,Ov 


+  Y;(e,e')  /  dTyLv(ty,e)exp[- 
o 


T  —  T  —  T 

Ov  V  V 

cose1 


]} 


and  the  spectral  sheet  emissivit  becomes: 


Tov  */2 

eSv  =  2e0  /  dTu/  d9  Si 
0  0 


3  Lv(fv,e)exp[- 


TOV-TV 

cose 


]  • 


(17) 


Because  of  our  assumption  that  t  »  optical  properties  of  the  droplet  medium 
are  frequency-independent,  the  si  ctral  variation  of  droplet  sheet  emissi- 
vity  arises  solely  from  the  freqi  ncy  (or  z)  dependence  of  $  from  which  y+ 
and  y~  are  computed. 

In  lieu  of  finding  solutions  of  eq.  (16)  using  Mie  scattering  phase 
functions  to  determine  yt ,  a  simi  ier  approach  was  used  to  find  the  magnitude 
of  the  effects  of  non-isotropic  scattering  on  e$.  We  neglect  the  frequency 
dependence  of  4  (or  equivalently,  simply  solve  eq.  (24)  at  a  single 
frequency)  and  represent  4  as  a  series  of  Legendre  polynomials: 

N 

4(a,o')  =1+1  P£(n*n')a£ 

4=1 

where  the  argument  is  the  cosine  of  the  angle  n  between  incident  and 
scattered  light,  and  the  leading  term  of  unity  (aQ=l)  is  required  to  insure 
that  /da  4(a,a')  =  4tt.  The  coefficients  a£  are  chosen  to  represent  various 


-15- 


types  of  scattering  (isotropic,  predominantly  forward,  predominantly 


backward,  and  perpendicular),  and  the  droplet  sheet  emissivities 


corresponding  to  these  scattering  types  are  determined.  This  representation 
greatly  reducs  the  computational  complexity  since  the  functions  y4  can  be 


evaluated  analytically: 


y±(9,0‘)  =  2tt  l  apP.(cose)Pf (±cosq ') . 


In  addition  the  intensity  function  L  can  be  represented  as  a  series  of 


Legendre  polynomials: 


L(t,u)  =  Lq(t)  +  l  L£(  f  )Pfi.(^  )a£ 
1 


where  y  =  cose. 


The  components  L»  obey  the  relations: 


N  1  H 

0)  r  r  dll 


l  (  /^VOexpt-IlL.  ] 

m=0  o  o 


+  /  dt  ‘ Lm( t ' )exp[- 


T0-T-T 


L£(t)  =  2™  l  f  P£(u)Pm(u)  {  /di ' Lm(T  ■  )exp[- 


m=0  o  p 


M-l )*/”  ]| 

0  M 


and  the  droplet  sheet  emissivity  is: 
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es  =  2eo  /  {Lo(T)El(  to"t  ) 
o 

N  1  T  _T 

+  l  a£Ls(T)  /^M^)  exP[-  — -  ]}  • 

£=  1  0  W 

The  integrals  over  u  can  be  evaluated  immediately  to  yield  a  series  of 
exponential  integrals,  since  the  functions  P£(u)  are  simple  polynomials  of 
degree  £.  One  is  left  with  a  series  of  integral  equations  which  are  solved 
by  discretizing  the  variable  t  and  solving  the  resulting  set  of  linear 
equations  by  matrix  inversion.  The  solutions  are  parameterized  by  the 
optical  depth  tq,  the  droplet  emittance  eQ,  and  the  coefficients  a^,...,aN. 

This  solution  procedure  was  carried  out  for  values  of  eQ  ranging  from 
0.1  to  1  and  for  optical  depths  ranging  up  to  4.0,  with  the  following 


choices  of  coefficients: 

a£>l  =  0 

isotropic  scattering 

al  =  a£>l  =  0 

forward  scattering 

al  =  _1»  a£>l  =  0 

backward  scattering 

a2  =  _1»  a£*2  =  0 

90°  scattering 

The  isotropic  scattering  case  corresponds  to  our  previous  results  (Ref.  1). 
Figure  4  shows  the  dependence  of  droplet  sheet  emissivity  on  tq  and  eQ  for 
isotropic,  forward  and  backward  scattering.  Forward  scattering  is  seen  to 
increase  e$  while  backward  scattering  decreases  es,  with  respect  to 
isotropic  scattering.  Emissivity  values  for  90°  scattering  were  found  to 
lie  between  the  values  for  forward  and  backward  scattering.  For  droplet 
sheets  with  the  product  e0t0  =  0.5,  non-isotropic  scattering  has  a 
significant  effect  on  sheet  emissivity  only  at  low  droplet  emittances  eq. 

An  emittance  eQ  =  0.1  and  an  optical  dpeth  t=5  yields  e$  = 


0.50  for  forward 
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scattering,  e$  =  0.44  for  backward  scattring,  and  es  =  0.46  for  isotropic 
scattering.  These  calculations  were  only  carried  out  to  determine  the 
maximum  effects  non-isotropic  scattering  might  have  on  sheet  emissivity.  In 
fact,  the  variation  of  the  angular  scattering  distribution  with  frequency 
makes  the  effect  of  non-isotropic  scattering  on  sheet  emissivity  much  less 
than  for  this  “worst  case"  model. 
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Fig. 


1  Baseline  configuration  for  the  liquid 
droplet  radiator. 


Fig.  2 


Relative  temperature  profile  in  a  droplet 
sheet  which  has  cooled  to  a  mean  temperature 
of  0.8  of  the  initial  uniform  temperature. 


optical  depth 


Droplet  sheet  emittance  vs  optical  depth  for 
non-isotropic  scattering  (Legendre  series  approx.) 


